Trimer States, Scale-invariance, and Effective Dimensionality 
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The momentum distribution is a powerful probe of strongly-interacting systems that are expected 
to display universal behavior. This is contained in the contact parameters which relate few- and 
many-body properties. Here we consider a Bose gas in two dimensions and show that, in spite 
of the presence of two three-body bound states, the two-body contact parameter is universal. We 
also demonstrate that the momentum distribution at next-to-leading order is vastly different from 
its three-dimensional counterpart. Based on this, we propose a scheme for measuring the effective 
dimensionality of a quantum many-body system by exploiting the functional form of the momentum 
distribution. 
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Introduction. Strongly interacting quantum systems 
are ubiquitous in Nature and naturally at the forefront 
of physics research. However, the theoretical study of 
strong interactions can be very difficult since our usual 
and intuitively clear perturbativc methods can fail mis- 
erably when inter-particle forces are strong. A success of 
cold atomic gas physics is the ability to create and manip- 
ulate strongly interacting gases in experiment One 
particularly interesting aspect is the ability to change 
the dimensionality by applied optical fields. This means 
that low-dimensional dynamics which is typically found 
in condensed-matter systems of great technological inter- 
est can be addressed. The held thus provides a testbed 
for models of strongly-coupled dynamics that are used to 
describe interesting materials. 

A breakthrough in the study of strongly-interacting 
quantum systems with short-range interactions was the 
derivation of a set of universal relations that relate 
the two-body correlations to the many-body thermody- 
namics through the so-called contact parameter, Ci [2|. 
These relations were subsequently confirmed in experi- 
ments on two-component Fermi gases [H, 0| and hold for 
bosonic gases as well [5|48| as confirmed by recent ex- 
periments Q. However, in the bosonic case, three-body 
correlations are very important and this implies that one 
must consider also a three-body contact parameter, C3 
0,1- 

A second avenue that is enjoying great success at the 
moment, is the experimental study of two-dimensional 
(2D) atomic Fermi gases Universal contact re- 

lations should also hold in this case SEMI (Ref. [3 
reports on recent experimental data concerning these is- 
sues). Interestingly, a recent experiment has found 
that the monopole breathing mode is essentially un- 
damped and has no interaction-dependent shift [20|-|22j. 
This implies a scale-invariance in the system [23| that has 
also been observed in weakly-interacting 2D Bose gases 



241 ], and which is hard to reconcile with the fact that a 
scale is provided by the energy of the two-body dimer 
which is always bound for attractive short-range interac- 
tion in 2D. 

From a few-body perspective, the scale-invariance of 
the 2D system is manifested in the spectrum of three 
identical bosons with attractive zero-range interactions. 
Here one finds that there are exactly two bound states 
which have energies E 3 = 16.52£ , 2 and E 3 = 1. 270^2 
with Ei the dimer energy [25l - |28| . This is in sharp con- 
trast to 3D, where an infinite set of geometrically sepa- 
rated states appear at the threshold for two-body binding 
(29| . In a real system, this scaling is broken by the finite- 
range of the interaction ^30] which ushers in the famous 
three-body parameter, ft 3 , . In 2D, scale-invariance is 
preserved at the three-body level, i.e. there is no k 2D 
and the energies must be proportional to the dimer en- 
ergy. 

In this paper, we study identical bosons in 2D with at- 
tractive short-range interactions and use few-body meth- 
ods to determine C2 and C3. This is achieved by com- 
puting the momentum distribution for three identical 
bosons, in particular its asymptotic behavior for large 
momenta. We provide both analytical and numerical ev- 
idence that support a universal tail behavior that is the 
same for both ground and excited states. This is the 
first time that C3 has been discussed in 2D to the best of 
our knowledge. Moreover, we show that the sub-leading 
term has a novel behavior that is radically different in 
2D as compared to 3D. Based on this fact, we propose 
to use the momentum distribution to measure the effec- 
tive dimensionality of a quantum system in the universal 
regime. Our study is thus a first step in exploring effects 
of dimensional crossover on higher-order correlations in 
many-body systems. 

Method. We consider three identical bosons with 
mass m. We use attractive two-body interactions of 
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zero range and parameterized by the dimer binding en- 
ergy, E 2 . The two-body T-matrix for energy E is thus 
t(E) = (— 2TT\n^/—E/E 2 )~ 1 in units where h = m = 1 
HI 31 , 32 1 • By using Faddeev decomposition and bosonic 
symmetry, the three-body wave function, VP, can be writ- 
ten 



, (q , p) _/{«>+/(|p-ll)+/(|p+ll) 



E 3 



p 2 + §q 2 



(1) 



where p = \{k\ — k 2 ) and q = |fea — + ^2) are 
Jacobi momenta, ki , i — 1,2,3 the lab momenta, and 
£3 is the three-body energy. The spectator functions, 
/(<?), satisfy the set of integral equations 



/(q) = 2r (i?3-|q 2 ) /d 2 * 



/(k) 



£ 3 



k 2 -kq' 



(2) 



Armed with the solution to this equation, the momentum 
distribution is 



n(q) = / dp 



/(<?)+/(|p-f|)+/(|p+f| 



ES 



p 2 + |q 2 



(3) 



Following the discussion in Ref. [8|, we define four com- 
ponents n m (q) = X)i=i n r(9)' where m = denotes the 
ground state and m = 1 the excited state. The individual 
components are 



<(<?) = f 2 m (v) / rf 2 P 
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n 2 n (q)=4f m (q) / rf 2 /c 



p 2 + |q 2 ) 



*7m(g) 

s 3 m + fq 2 : 



Z' TPT. 



n™(q) = 2 d 2 k 



k 2 + q 2 



kqr 



= 2 / d 2 fc 



(£™ + k 2 + q 2 + k • q) 



/m(fc)/m (|k - q|) 

(£" + k 2 + q 2 + k ■ q) 



(4) 
(5) 

(6) 

(7) 



where m on /(g) and E 3 labels the state. Throughout, 
we measure all momenta in units of \JE 2 . 

Large-momentum limit. The leading order (LO) be- 
havior of the momentum distribution exhibits the same 
C 2 k~ A tail in ID, 2D, and 3D since it derives solely from 
two-body physics [33(. However, C2 depends on what 
system is addressed and whether few-bodybound states 
are present. For bosons in 3D, the tail is 0, Q 



1 cos[25 O ln(V3fc/K,) + 0] ^ . . 
n 3 D{k) -> -^62 H p C 3 , (8) 



where so = 1.00624 and <fi — —0.87280 are constants that 
can be determined from a full solution of the three-bosons 



problem in 3D at unitarity @ with trimer energy E 3 = k 2 
(using = for simplicity). The log-periodic three- 
body next-to-leading order (NLO) term derives from the 
Efimov effect, whose solution can be used to determine 
C 2 = 53.097/ft* and C 3 = -89.263/ft 2 [|. As discussed 
above, in 2D there is no Efimov effect for three bosons. 
The log-periodic behavior is therefore not expected a pri- 
ory. As we will now demonstrate, the distribution in 2D 
is very different. It has the structure 



n 2D {k) -> -£jC 2 



ln 3 (fc) 



(9) 



and we see indeed a very different NLO term. This 
implies that the NLO term is in fact an effective mea- 
sure of dimensionality of bosonic systems in the universal 
regime. We will return to this point below. 

To derive the tail behavior in Eq. ©, one needs to 
determine first the spectator function, f m {q), in Eq. ([2]) 
for large q. This can be done analytically and we pro- 
vide the details in the supplementary material [34| . The 
result is that f m {q) A m \a(q)/q 2 , where A m is a state- 
dependent constant. This function can now be inserted 
into Eqs. (J4j) to {[TI) and the momentum tail can be de- 
termined. The technical details are given in (34|. After 
the dust settles, the tail behaviors can be written 



n?(q) 



4irA 2 m ]n 2 (q) 



n 2 n (q) -> 4^ 



A 2 „M 3 (q) 



— dkkf 2 m {k),n^{q) 
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A 2 n \n 3 (q) 



The LO term clearly comes from n™(q), while NLO has 
contributions from ri 2 n (q) and n™{q). However, there is 
an additional complication as NLO will also come from 
rig 71 at the next order (not shown above) . More precisely, 
we need to determine 



D r , 



= lim 

q— >oo 



ln 3 (gV 



(10) 



where C 2 = lim^oo q n m {q) which is independent of m 
as we discuss below. We find that D m is a non-zero con- 
stant that depends on the state m, which means that C3 
in 2D should be denoted C™. Adding the NLO contri- 
butions from 
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1™ (q) , n™{q) and n™(q), we find 
1.01 



C° 3 



52.07 and Cg 1 



(11) 



This m-dependence is absent in 3D for k* — > but at 
the cost of the log-periodic term due to the Efimov effect 
[8j. That the present 2D case has state-dependence is a 
result of the lack of geometric scaling symmetry in 2D. 
Note that the next order comes from n™(g) and differs 
by one power of ln(g) compared to the NLO term. 

Universal behavior. The LO behavior in 2D is char- 
acterized by C 2 . Explicitly, we have 



»§(?) 



3.71E 2 



and n\(q) 



0.28E 2 



(12) 
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FIG. 1: LO momentum distribution tail, q 4 n(q), for ground 
(upper solid black line) and excited (lower dashed red line) 
three-body states. Note that the vertical axis is not uniform. 
The asymptotic dashed lines are the analytical results [34|. 



These results have been obtained analytically (see [34|). 
We have also done a numerical check which is shown in 
Fig. [I] The units in Eq. (|12[) are, however, not natural in 
the same way that is seen in Eq. ([5]) where k» provides 
the overall scale. The natural scale is E3, and using this 
we find 3.71£ , 2 /16.52£; 2 = 0.224 and O.28E2/I.27OS2 = 
0.219 for ground and excited state respectively. This is a 
striking results that demonstrates the state-independence 
of the LO term in 2D to within our numerical accuracy 
of about 2%. We thus predict that the two-body contact 
for a bosonic system in 2D with short-range attractive 
interactions is 



C 2 /E 3 = 0.222 ± 0.0025, 



(13) 



where E3 is the trimer energy. 

The universal tail behavior is far from trivial. In 3D 
and at unitarity, the discrete scale invariance induced by 
the divergence of the three-body problem, implies that 
the system should behave similarly irrespective of which 
trimer state one considers. This does not occur in 2D 
and the universal trimer energies are in some sense magic 
numbers multiplying the only scale available, E 2 - Our 
results show that in spite of this major difference, the 2D 
momentum tail displays universal behavior, i.e. C2/E3 
has the same value for both ground and excited states. 

Dimensional crossover. Comparing the expressions 
in Eq. ([5]) and Eq. ([5]), we see the same LO behavior 
at large momenta, but a vastly different NLO term. The 
oscillations seen in Eq. (|8]) can be traced directly to the 
discrete scaling symmetry, or more precisely, the break- 
down of scale-invariance in the system. It is known that 
the condition on the dimension, D, for this behavior is 
2.3 < D < 3.8 27J, |28[. If we imagine an interpolation 



between 2D and 3D, we would expect to see log-periodic 
terms in this range of D. The NLO term is therefore a 
tell-tale sign of effective dimensionality of the system as 
we will now discuss. 

In experiments that study cold 2D quantum gases, one 
uses a tight transverse optical lattice potential to reduce 
the motion in this direction [1]. As recent experiments 
have beautifully demonstrated, the strength of the trans- 
verse optical lattice can be used to interpolate between 
2D and 3D behavior of fermionic two-component systems 
1 lj . Here we are concerned with bosonic systems, 
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and our results above demonstrate how one can use the 
tail and in particular the NLO part of the momentum dis- 
tribution as a measure of the effective dimensionality felt 
by the particles in the system by identifying the presence 
of log-periodic behavior. In Fig. [2] we show the extreme 
cases of 2D and 3D where the log-periodic oscillations are 
clearly seen in the latter, while the former has a smooth 
behavior. 




FIG. 2: NLO momentum distribution comparison of 3D (up- 
per panel) and 2D (lower panel). The 2D momentum distri- 
bution is the one of the ground state, but the result is similar 
for the excited state. 



A measurement of the overall functional form of the 
NLO term is thus enough to determine the effective di- 
mensionality of the squeezed bosonic gas. In a real ex- 
periment, the motion in the transverse dimension is of 
course quantized by the lattice, and to get a full quanti- 
tative understanding this must be taken into account (see 
for instance Ref. [35j). However, since experiments have 
shown that it is possible to reach both the extreme 2D 
and the 3D regime, there must necessarily be a dimen- 
sional crossover that can be seen in the NLO behavior. Of 
course, from a theoretical point of view it would be very 
attractive to be able to map the strength of the trans- 
verse confinement into some effective dimensionality D e g 
which could be non-integer 33]. 
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Experimental implementation. As we have demon- 
strated, the NLO term in the momentum distribution 
carries a tell-tale signature of the dimensionality of the 
quantum system under study. The 2D to 3D crossover 
is of immense interest at the moment [ID, (3||. To 
connect the results obtained here to actual experiments, 
we need to consider our units, the dimer binding energy 
£2, and the effects of the transverse confinement on this 
two-body bound state. 

The interaction is controlled by Feshbach res- 
onances (3(|. However, under the confinement, 
the dimer energy is modified and becomes E2 = 
Bftijj z exp(—y/27rl z /\a3B)/'K 37]. Here lu z is the trans- 
verse harmonic confinement frequency, l z = y/ h/mui z the 
trapping length, a the 3D scattering length associated 
with the Feshbach resonance, and B — 0.905 is a con- 
stant. This formula holds for a < and \a\ <C l z , while 
on resonance, \a\ — > 00, E2 — Q.2AATujj z . Corrections arise 
from the non-harmonic optical lattice [38| . but they are 
not essential for our purposes. The dimer energy scale 
can be converted into a momentum scale, through 
E2 = h 2 k^/2m. To access the tail behavior and the 2D- 
3D crossover, we see from Figs. [T] and [2] that the range 
k ~ 10 1 — 10 3 fco is sufficient. Recent 2D Bose gas exper- 
iments 
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|39|, use l z ~ 3800ao, where ao is the Bohr 
radius, which implies that ko ~ 10 _4 Oq 1 when la]_ = 00. 
For the momentum distribution measurements [3j, |4| , the 
maximum momentum reported is about k ~ 10~ 3 a^~ . 
This implies that an order of magnitude or two beyond 
the reported capabilities is necessary. However, if a is 
tuned away from resonance to the a < side, E2 will de- 
crease rapidly according to the formulas above, inducing 
a corresponding rapid decrease of fco which should render 
the physics discussed here within reach of current experi- 
mental setups. Notice that the van der Waals length scale 
of about lOOao is in the deep tail, so there is no conflict 
with the universal zero-range description employed here. 

Summary and Outlook. We have taken a first step in 
the study of higher-order correlations and dimensional 
crossover by demonstrating how trimer observables in 
strongly-interacting quantum gases can be used to probe 
dimensionality. Specifically, we see the breakdown of 
scale-invariance directly in the functional form of the tail 
of the momentum distribution. 

A clear direction for future study is a full inclusion 
of the transverse direction and the discrete spectrum it 
brings. We have shown that a crossover with fundamen- 
tal influence on the momentum tail will happen, but map- 
ping it out in a system that is squeezed by optical lat- 
tice^) is the next task. This would also be interesting 
for the 1D-3D or 1D-2D crossovers. Another venue to 
explore is mass imbalanced systems where the spectrum 
is known to be more rich in 2D than the equal mass case 
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SUPPLEMENTARY MATERIAL 
MOMENTUM DENSITY ASYMPTOTIC BEHAVIOR 

In this supplementary material we give the technical details of the analytical and numerical determination of the 
three-body wave function and momentum distributions. We will use units h = m = 1 and all energies are given in 
units of the two-body dimcr energy E 2 implying that all momenta are in units of \/Ea. The one-body momentum 
distribution is given by 



n 



/(9)+/(ip-fi)+/(ip+!r 1 



E 3 + P 2 + |q 2 



47T „ , N , / / 3 „ \ 47T 



(14) 



i(g) = J d 2 p 
and can be split into four parts through 

4 

n m {q) = Y.<^ ( 15 ) 
1=1 

where the subscript to distinguish ground (to = 0) and excited states (to > 0). The individual expressions are 

n?( q ) = f\ q ) I d 2 p — 1 = J^L, (16) 

J (E 3 + p 2 + f q 2 ) H + 4^ 

<(*) = m)I#P { j}Xf^ 

n?(<l) = 2 / d V /2(|P 2 + !'L = 2 / d-k -J^- .„ (18) 
J (£ ; 3 + P 2 + |q 2 ) J (£ 3 + k 2 + q 2 +k-q) 

n m (a) - 2 f j2 /(|p+f|)/(|p-fl) 

= 2/-d 2 fc , F / it ) /i |k 2 " q J ) , a • (19) 
Here we are interested in the limit q — > oo, where we find 

»T(«)»Y^. (20) 



<(g)«^/ 2 ( < z)ln h/i?3 + i q 2 I «^/ 2 («)ln(?). (21) 



4 /■oo 

« -4 / *fc/ 2 (fc), (22) 
1 Jo 



The asymptotic form of /(g) when q — > oo is (see derivation in the next section) 

/(g) ^(j-s-oo A m ~~^2~- 



(23) 



(24) 



The functional form of asymptotic behavior is the same for both ground and excited states. However, the normalization 
constant, 



A m = lim f m (q) — , 
g->oo in q 

is different. Using the condition J d 2 k |n(fc)| 2 = 1, one finds 

A Q m 0.236 and A 1 w -0.103 
The numerical results for the asymptotics are shown in Fig. [3] 



(25) 



(26) 





FIG. 3: Asymptotic behavior of f(q) for the ground state (left) and the first excited state (right). 
Inserting Eq. (j!M)) ) in Eqs. (|20l) )- (|2"3")) one obtains the normalized asymptotic behavior 



, Air ln 2 (g) 
3 g 6 ' 



<(g) 



with C m = / dkkf(k), 

g /o 



(27) 
(28) 
(29) 

(30) 



The normalization constants given in Eq. (I26[) determine the asymptotic values and behaviors of the partial mo- 
mentum density in Eqs. ()27|) - (I30I) when q — > oo which are 



i 2/ \ 

n?(g) ->• 13.50 — and nj(g) -> 0.26 



ln 2 (g) 



(31) 
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n° 2 (q) -> 40.50 ancl n \{q) -> 0.78 



(32) 



o, 3.71 x , 0.28 



(33) 



n2(g) -> 20.25 and n J( 9 ) ^ .39 



(34) 



ASYMPTOTIC FORM OF f(q) 

In the symmetric case, where the three masses and the three two-body binding energies are set equal to one, the 
spectator function fulfills the integral equation 



/(?) 



1/tt 



d l k- 



f(k) 



In U Ez + 3/AqA J E 3 + q 2 + k 2 + k-q 



(35) 



This can be cast into the useful form 



/(<?) 



1/tt 



dk 



kf(k) 



In (y/E 3 +3/4q 2 ) Jo E 3 + q 2 + k 2 J Q 1 + a cos 6 ' 
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d0 



(36) 



with a — kqj (E3 + q 2 + k 2 ). The angular integral is 



de 



2tt 



and one obtains 



/(g) 



o 1 + a cos 9 y/l - a 2 



dk 



for < a < 1 



kf(k) 



In (V^ + W) Jo (Es + q 2 + k 2 )^l- 



which can be rewritten as 



(37) 
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I A dk - 



kf(k) 



(E 3 +q 2 + k 2 )Jl 



fc/(fc) 
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(39) 



where A is a large-momentum cut-off that will be useful below. Taking A >> VE3, the spectator function in Eq. (|39 
is approximately given by 



/(g) 



Jo A dk ' 



In(VBs+3/V) J0 (E3+s2+fe2) XT 



/a 00 *- 



fc/(fc) 



ln(V-E3+3/V) JA (g 2 +fe2) 



(40) 



For q — > 00, the first term on the right-hand-side of Eq. (1401) tends to zero in the following manner 



q 2 In q 



*/(*) 



(41) 



8 



Now, we assume that 



/(?) 



In q 



(42) 



Inserting this ansatz and taking the limit q —¥ oo in the second term on the right-hand-side of Eq. (|40|) one finds 



Ha) 



In q 



dk 



k\nk 



dk-. 



Ink 



when q — > oo. Changing variables to y = k/q, the second spectator function term becomes 



h{q) 



dy hxy + hxq 
q 2 \nqJ A/q y (1 + y 2 ) 



which can be split into 



Mo) 



q 2 In q 



dy lay 

A/ 9 7(i + y 2 ) 



In q 



dy 



a/, i/ (i + y 2 ) 



The first integral term on the right-hand-side of Eq. (|451) is 



ciy lny 



1 In 2 y 



A/q 



/a/, 1/(1 + ?/ 2 ) 2 1 + y 2 
The second term on the right-hand-side of Eq. (|4"5")) is 

dy 1 In y 



yln y 



A/ q (i + y 2 ) 



2\2 _ 1^°° 



■In 2 



A 



■ In 2 g. 



A/ 



, y (l + y 2 ) 1 + y 2 



A/q 



dy 



y\ny 



A/q "(i + y 2 ) 



2\2 ~ 9- >oc 



In 



A 



lng 



Inserting the results of Eqs. (|4l))) ) and (|47l) in Eq. (|4^)l ) we arrive at 

/a (9) 



(In g- - In q) = —. 



q \nq 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



Collecting the results Eqs. (|4T|) and (|4^|) we conclude that the ansatz in Eq. (l4"2"j) give us the asymptotic behavior of 
the exact spectator function. 

In Fig. Uit is showed both the spectator function obtained from the numerical solution of the set of coupled integral 
equations and the spectator function asymptotic behavior given in Eq. (|42p for the ground state. The log- log scale show 
us that both magnitude and line's inclination are very close in the region of q between 100 and 2000 for the numerical 
and analytical calculations. For q > 2000, the curve which represents the numerical solution of the integral equations 
starts to oscillate around the analytical form. 
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FIG. 4: Spectator function, /(g), for the ground state calculated numerically (black solid line) and using the ansatz f{q) = Ao^t' 
(red dashed line). The solid (black) line tends to oscillates around the dashed (red) one as q — > oo due to finite numerical 
precision. 



